Group theory analysis of electrons and phonons in N-layer graphene systems 
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In this work we study the symmetry properties of electrons and phonons in graphene systems as 
function of the number of layers. We derive the selection rules for the electron-radiation and for the 
electron-phonon interactions at all points in the Brillouin zone. By considering these selection rules, 
we address the double resonance Raman scattering process. The monolayer and bilayer graphene in 
the presence of an applied electric field are also discussed. 
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I. INTRODUCTION 

The current interest on graphene and its multilay- 
ered materials has been stimulated by various experi- 
mental and theoretical works addressing the physics of 
Dirac fermions and the potential for device apphcations 
( [H) 0| and references therein) . Group theory is a power- 
ful theoretical tool to determine eigenvectors, the num- 
ber and the degeneracies of eigenvalues and to obtain 
and understand the selection rules governing, for exem- 
ple, electron-radiation and electron-phonon interactions. 
Although the symmetry aspects of mono-layer graphene 
and graphite have been largely discussed in the literature 
, the recent findings generate interest in a group theory 
analysis depending on the number of graphene layers. 

This work presents a group theory analysis for elec- 
trons and phonons in mono-, bi- and tri-layer graphene, 
extending for N layers depending if N is even or odd. The 
selection rules for electron-radiation interaction within 
the dipole approximation and for electron scattering by 
phonons are derived. With these selection rules, we dis- 
cuss the double-resonance Raman (DRR) scattering pro- 
cess, which has been widely used to characterize the num- 
ber of layers d, H, Q and toprobe their electronic and 
vibrational properties 0, H, [^. Finally, we also discuss 
the differences when mono- and bi-layer graphene are ex- 
posed to external electrical fields, giving insight on the 
gap opening in the biased bilayer graphene [Tol. Till. fl^. [Tsj 
and different selection rules for the electron-phonon scat- 
tering (EPS) process. 

Section |TT] gives the symmetry properties for mono- 
layer and N layer graphene depending if N is even or 
odd. The monolayer and bilayer graphene are also con- 
sidered in the presence of an electric field perpendicular 
to the graphene plane. The notation adopted is related to 
the space group symmetry, and conversion to the point 
group notation can be found in the appendiji, Section 
mil presents the selection rules for the electron-radiation 
interaction. Section ITVl shows the V point Raman and in- 
frared active modes, and we extend the electron-phonon 
selection rules to points in the interior of the Brillouin 
zone in section |Vl Considering both sections IIIII and |V] 
we address the DRR process for mono-, bi- and tri-layer 
graphene in section IVIl The main findings are summa- 



rized in section IVTIl 



II. SYMMETRY PROPERTIES 
A. Group of wavevector 

Figure [Ha,d) shows the hexagonal real space for the 
monolayer graphene with two inequivalent atoms in the 
unit cell. The origin is set at the highest symmetry 
point, i.e. at the center of a hexagon. The recipro- 
cal space is shown in Fig. [IJg) highlighting the high 
symmetry points F, K, K', M and lines T, T', E. Any 
other generic point outside the high symmetry lines and 
points is named here u. The monolayer graphene on an 
isotropic medium has the space group P6/mmm (Dg,J 
in the Hermann-Mauguin notation. At the F point, the 
group of wavevector (GWV) is isomorphic to the point 
group Dgh (the Schoenflies character tables for the point 
groups can be found in Ref. [1|). 

The real space for bilayer and trilayer graphene with 
AB Bernal stacking are show in Figs. [Hb,e) and (c,f), re- 
spectively. The symmetries for N-layer graphene, with N 
even or odd (from now on, 1), are the same of bilayer 
and trilayer graphene, respectively. The main symmetry 
operation distinguishing the point groups between even 
and odd layers are the horizontal mirror plane, which 
is absent for N even, and the inversion, which is absent 
for N odd. The point groups isomorphic to the GWV 
for mono-, N-layer graphene (N even and odd), and for 
N infinite (graphite) are listed in table[T] for all points 
and fines in the first Brillouin Zone (BZ). The GWV for 
N-layers graphene are subgroups of the GWV for single 
layer graphene. The direct product between the group 
from N even and N odd gives the graphene GWV, i.e. 



{Ge.en|0}®{God<i|0} = {G 



monolayer 



|0}. 



On graphite, the wavevector point groups are isomorphic 
to the wavevector point groups of monolayer graphene, 
but differ fundamentally for some classes where a transla- 
tion of c/2 is present, graphite belonging to the mmc 
(Dg^) non-symorphic space group. 
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FIG. 1: (a) The real space top-view of a monolayer graphene showing the non-equivalent A and B atoms, (b) The real space 
top- view of bilayer graphene. Light and dark gray dots, and black circles and dots represent the atoms in the upper and lower 
layers, respectively, (c) The real top-view trilayer graphene. Small and large black dots, light and dark gray dots, and traced 
and solid circles represent the atoms in the upper, middle and lower layers, respectively. The unit cell of the (d) monolayer, 
(e) bilayer and (f) trilayer graphene. (g) The reciprocal space showing the Brillouin zone in light gray, the high symmetry 
points and lines and the two primitive vectors. 



TABLE L The space groups and wavevector point groups for 
mono-, N-layer graphene and graphite at all points in the BZ. 





Space group 


r 


K ,K-. 


M 


T (T') 


s 


u 


Monolayer 


P6/mmm 




D3h 


Dah 


C2u 


C2i> 


Cih 


N even 


P3ml 


D3d 


D3 


C2h 


C2 


Cii, 


Ci 


N odd 


P6m2 


D3h 


Csh 


C211 


Cih 


C2i> 


Cih 


N infinite 


P63 jmrac 




D3h 


D2h 


C2u 


C2d 





B. Lattice vibrations and tt electrons 

The representations for the lattice vibration (Fiat. nib.) 
and for the tt electrons (Ftt) are given by ^lat.vib. = 
pe? (g, p^ector ^ pe, ^ p^ ^ respectively, where 

r*^' is the atom equivalence representation, p'"^'^*"'' is the 
representation for the vectors y and z. For F^ we used 
only F^, which is the irreducible representation for the 
vector z, since tt electrons in graphene are formed by p^ 
electronic orbitals. The results for all points and lines in 
the first BZ for the Viat.vih. are found in TablelTTl and for 
the F^ in TableHIIl 

Table IIIII shows that the tt electrons in monolayer 
graphene are degenerated at the K (Dirac) point, as ob- 
tained by theory [l^ . Figures [21^ a) , (b) and (c) show 
the electronic structure of a mono-, bi- and tri-layer 
graphene, respectively calculated via density functional 
theory (DFT) [l^ [3, [13] ■ The symmetry assignments 
of the different electronic branches shown in Fig. [2] were 
made according to the DFT projected density of states. 

The bilayer graphene with AB Bernal stacking [see 
Fig.[IJb,e)] is also a zero gap semiconductor composed 
by two conduction and two valence bands, and the elec- 
trons exhibit a parabolic dispersion near the K point. 



Two bands are degenerated at the K point (see Table|TTT] 
and Fig. [Ifb)) and the other two have a gap of 271 , where 
71 is the Slonczewski-Weiss-McClure parameter [isl Il9ll 
that have a experimental value of ^ 0.3 - 0.4eV [7l.llCll.[20l|. 

Trilayer graphene in the ABA Bernal stacking (see 
Fig. [Hc,f)), belongs to the D3/1 point group and Fig. 
[2l[c) shows its electronic dispersion. The K point of 
trilayer is isomorphic to C3/1. In Tables [Til and IIIII 
and K+* are the two one-dimensional representa- 
tions of the K'^ representation, where * means the com- 
plex conjugate. The same happens for the repre- 
sentation. The electron representations will be given 
by Ff = X+ -I- 2Xf -f K^* + K:^ -h K:^* for the K 
point and Ff' = A'+ + + X+ + K:^* + if 2" for K' 
point. Although time reversion symmetry can imply de- 
generency between complex conjugate representations in 
ciclic groups, in graphene the complex conjugation also 
takes K into K' point and, consequently, there are no de- 
generated bands at the K (K') point, in agreement with 
tight-binding calculations when including the 72 and 75 
next-nearest-layer coupling parameters [2l|, [g^] . This en- 
ergy gap is also obtained from ab initio calculations (see 
the inset of Fig. H (c) and Ref. [Hi). 



C. Gated mono- and bi-Iayer graphene 

If the monolayer graphene is in the presence of a per- 
pendicular electric field (gated graphene), the Fermi level 
changes. The presence of charge inhomogeneity caused 
by substrate and/or absorbed water can generate the 
same effect as verified in transport [U, [2^ and Raman 
measurements [26l| , where the Dirac point is shifted from 
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TABLE II: The Fiat.vib. wavevector point-group representations for mono- and N-layer graphene at all points in the BZ. 





Monolayer 


N even 


N odd 


r 


+ + r+ + r+ 


Ar(r+ + r+ + r- +r-) 


(AT - i)r+ + {N + i)r- + {N + i)r+ + (at - i)r- 


K 


Ki +K+ +K+ + 


N{Ki +K2 +2K3) 


+ + [/(JV) + 2]K+ -1- IfiN - 2)]K+-' 
+NK2 +(N -1)K2' " 


M 


M+ + M+ + M+ + M- + M- + M- 


Ar(2A/+ + M+ + M," + 2M.7) 


2NMi + {N - 1)M2 + (N + 1)M3 + 2AfM4 


T(T') 


2Ti +T2 + 2T3 + T4 


3Ar(Ti + T2) 


(3W + 1)10T+ + [ZN - 1)T" 


S 


2Si + 2S3 + 2S4 


Ar(4Si +2S2) 


2NY.I + {N - 1)S2 + {N + 1)S3 + 27VS4 


u 


4ji+ + 2u^ 


6Nu 


(ZN + l)u+ + (3W - l)u- 



"Where f(N) = Em=o[©(^ " 4m - 2) + 3e(Af - 4m - 4)], where 
0(2:) is equal to if x < and equal 1 otherwise. 



TABLE III: The wavevector point-group representations for mono- and N-layer graphene at all points in the BZ. 





Monolayer 


N even 


N odd 


r 


r-+r,+ 


N{rt + r-} 


(AT - i)r+ + (jv + i)r- 


K(K') 







i^)K+ + {^)K- + g{N)K+'(K+) + g(N - 2)(K+)(K + *) + g(N)K- + g(N + 2)K-' " 


M 


Mg^ + M2 


N{M^ + 


(N - l)Mi + (N + l)Mi 


T(T') 


T2 + Ti 


Af(Ti +T2) 


(N - 1)T+ + (N + 1)T" 


E 


2S4 


2AfEi 


(TV - l)Si + (Ar + 1)S4 


u 


2m" 


2AfM 


(JV - l)ti+ + {N + l)u- 



"Where g{N) = Em=o ©(^ - 4m - 2), where e{a;) is equal to 



if a; < and equal 1 otherwise. 



the neutrality point. In this case, the tt electrons loose 
the horizontal mirror plane and the inversion symmetry, 
and the system is isomorphic to the point group Cg^. 
The irreducible representations for the for the gated 
graphene can be found in Table llVl There is no gap open- 
ing at K point for a perfect perpendicular electric field 
effect. 

The biased bilayer graphene have attracted a lot of at- 
tention recently because it is the only material known to 
have a tunable energy gap [13, [HI [13, [13] , promising for 
applications on devices and lasers with tunable energy. 
The mechanism behind this feature is based on applying 
an electric field perpendicular to the graphene layers, so 
that the two layers will be under an inequivalent poten- 
tial. Then it is possible to open a gap at the K point, 
breaking the double degenerated K3 irreducible repre- 
sentation into two one-dimensional irreducible represen- 
tations. Since the biased bilayer graphene brakes the 
inversion center symmetry, the group of the wavevector 
at r for perfect perpendicular electric field is isomorphic 
to TablellVl shows that the biased bilayer contains 
the two one-dimensional representations K2 and Kj at 
the K point, then a gap opening is expected on the basis 
of symmetry arguments of inequivalent layers. 

The representations for the Tiat.vib. of the gated mono- 
layer (or biased bilayer) are the same as the monolayer 
(or bilayer) in an isotropic medium, given in table [Til 
The electric field does not affect the symmetries of the 
phonons. 



III. SELECTION RULES FOR 
ELECTRON-RADIATION INTERACTION 

The symmetry properties described in the previous sec- 
tion will now be applied to physical processes. In this sec- 



TABLE IV: The group of wavevector and its representa- 
tions for gated monolayer and biased bilayer graphene. 





Gated monolayer 


Biased Bilayer 




GWV 




GWV 




r 




Ti + r4 




4ri 


K(K') 




Ks 


C3 


2Ki +K2 + K* 


M 


C2v 


Ml + Ms 


Civ 


4Mi 


T(T') 


Cii, 


T1+T2 


Ci 


4T 


S 


Cit, 


2Si 




4Si 


u 


Ci 


2u 


Ci 


4m 



tion we discuss the selection rules for electron-radiation 
interaction in the dipole approximation, with emphasis 
on the high symmetry lines T and T' in the electronic 
dispersion, where interesting phenomena occur. 

In the dipole approximation, the absorption of light 
in a material is related to the wave functions of the elec- 
tron states in the valence {ip"{'k)) and conduction ('0'^(k)) 
bands and the polarization of the incoming light (P) by 
W{k) oc |P • (^A=(k)|V|V'"(k))|2 [13,111. Knowing the 
symmetry of the initial an final states, and the repre- 
sentation that generates the basis function of the light 
polarization vector (x, y or z), group theory can be used 
to compute wether iy(k) is null or not. The results are 
summarized in TablelVl considering graphene layers lay- 
ing in the {x,y) plane and light propagating along z. In 
the case of graphene, the light absorption up to 3 eV 
occurs only at T, T' and u points. 

It is important to highlight some results given in Ta- 
ble [V] In the case of monolayer graphene on an isotropic 
medium, numerical calculations show an anisotropy in 
the optical absorption [H, HI, [lO] • This anisotropy has 
indeed a symmetry basis, as clearly seen when analyzing 
the selection rules at the T line. Absorption by visi- 
ble light has to couple T2 and T4 tt electron symmetries 
(see Fig.[2l^a)). For the T line direction along y, the only 
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FIG. 2: The electronic dispersion for the tt electrons calculated by DFT and the irreducible representations (F^) for (a) 
monolayer, (b) bilayer and (c) trilayer aloiigthe KFMK directions. The calculation was done via pseudopotencial DFT [15j] as 
implemented in the SIESTA program p^. [l7(|. We used a basis set composed of pseudo atomic orbitals of finite range and the 
Local Density Approximation (LDA) with the Ceperley- Alder parametrizationfor the exchange-correlation functional. 



TABLE V: Selection rules for electron-radiation interaction 
with X and y light polarization in mono-, bi- and tri-layer 
(see Fig. [1] (g) for x and y definition) . For N even and N 
odd the selection rules are the same as for bi- and tri-layer 
graphene, respectively. 





BZ point 


polarization 


W(k) 


monolayer 


T 


1/ e Ti 


Tb ig) Ta (g) T4 non null 
T2 ® Ti (g) T4 null 


u 


X, y g 


u~ (gi 11+ Cgi non null 


gated 
monolayer 


T 


X e Ta 


Ti ig) r2 (g T2 non null 
Ti ® Ti ig T2 null 


u 


X, y a u 


ii (g) n ® M non null 


bilayer 

(N-even) 


T 


X €T2 


Ti ig r2 (g Ti null 

Ti (g r2 (g T2 non null 

T2 Ig T2 «) T2 null 

Ti Ig Ti Ig Ti non null 

Ti Ig Ti Ig T2 null 

72 Ig Ti Ig T2 non null 


u 


X, y £ u 


u (g ji ig u non null 


biased 
bilayer 


T 


X, y eT 


T®T ®T non null 


u 


X, y £ u 


u (g u ig M non null 


trilayer 

(N-odd) 


T 


x,y eT+ 


r+ Ig r+ Ig T+ non null 
r+ Ig r+ Ig T" null 
T" (g T+ (g T" non null 


u 


X, y € 


It"*" ig u'*' ig li^ non null 
ti"*" Ig tt"*" ig u" null 
«~ (g Ji+ (g non null 



allowed absorption is for light polarized along the x di- 
rection. For incident light polarization along the y direc- 
tion, no absorption will occur along KF direction, giv- 
ing rise to the optical absorption anisotropy on graphene 
[28, 29, 30]. Outside the high symmetry T line there is a 
non-zero probability of absorption and the anisotropy is 
obtained by defining orthogonal basis, as shown in Ref. 

When the monolayer graphene is on top of a substrate, 
with the influence of the environment changing the Fermi 
level, there will be no change in the selection rules for 
electron-radiation interaction. Along the T line, the tt 
electrons are described by Ti and T2 representations, 




FIG. 3: (a,b) Schematic electron dispersion of bi-layer 
graphene along the KF direction showing the possible transi- 
tion induced by (a) a photon with T2 symmetry [x polariza- 
tion) and (b) a Ti photon (y polarization), (c) The electronic 
dispersion of tri-layer graphene showing the five possible tran- 
sitions by light absorption. 



where T2 and Ti contain x and y basis function, respec- 
tively. Again there will be no absorption for y polariza- 
tion. 

The bilayer graphene is composed by four electronic 
bands at the T line, belonging to two Ti and two T2 ir- 
reducible representations. The four possible transitions 
are illustrated in Fig. [3] (a,b). In this case both x and y 
polarized light can be absorbed. For the biased bilayer 
graphene, all electronic representations are the same, and 
it contains both x, y base functions for light polarizations. 
Thus, all the four transitions are allowed connecting all 
the four bands by the same light polarization, differently 
from the unbiased bilayer case where the light polariza- 
tion selects the pair of bands that can be connected. 

The trilayer graphene will have more possibilities for 
light induced transitions, since there are more possibil- 



5 



ities between the three tt and three tt* bands. Along 
T(T') direction, there are two T+ and four T~ bands 
giving raise to five possible transitions (see Table |V| , as 
shown in Fig. [3] (c). 



in the process. The allowed electron-phonon scattering 
processes for monolayer, gated monolayer, bilayer, biased 
bilayer an trilayer graphene along the KF and KM direc- 
tions (T and T', lines respectively) and at a generic u 
point are summarized in Table IVII 



IV. SELECTION RULES FOR THE 
FIRST-ORDER RAMAN SCATTERING AND 
INFRARED ABSORPTION PROCESSES 

The first-order Raman scattering process is limited to 
phonons at the center of BZ (F point) due to momentum 
conservation requirement (phonon wavevector q — 0). 
In monolayer graphene the first-order Raman spectra is 
composed by the G band vibrational mode, which is dou- 
bly degenerated at the F point with Tg symmetry. The 
Raman active modes depending on N (N>1) (without 
acoustic modes) are: 



^Raman 



= 7V(F+ +F+), for N even 



Raman 



Nr+ + {N- 1){T^ + F+), for N odd 

For even number of layers the G band belongs to the 
F;^ . There is a low frequency Fj^ mode with frequency de- 
pending on the number of layers (35-53 cm~^)[3l|. Two 
new Raman active modes near ~80 cm~^ and ~ 900 
cm^^ appear with irreducible representations [3ll.[3^. 
For odd number of layers the G band is assigned as a com- 
bination of and F;^ representations, also the smaller 
wavenumber component is active in Raman by a F^ rep- 
resentation. 

For monolayer graphene there is an infrared (IR) active 
mode belonging to the F^ representation, giving rise to 
an absorption near ~ 870 cm~^. The IR active modes 
for N>1 are: 



F^^ ^{N - l)r- + {N - 1)T-, for N even 
F^^ ^ N{T+ + (F2"), for N odd 

For even number of layers the active modes belong to 
the F^ and F;^ representations, the later one referring 
to the ~ 1590 cm~^ frequency vibration [3l|, |33|. The 
infrared active modes for odd layer number belong to T2 
and Fj^ which are also Raman active. 



V. ELECTRON SCATTERING BY g / 
PHONONS 

The electron-phonon scattering (EPS) is calculated 
from the initial and final electron wave functions coupled 
by the phonon eigenvector [H, Q using the phonon- 
induced deformation potencial. Therefore, the selection 
rules of the EPS processes are obtained by the direct 
product of the symmetries of the initial and final elec- 
tronic states and the symmetry of the phonon involved 



TABLE VI: Allowed processes for electron-phonon scattering 
for mono-, hi- and tri-layer graphene along the T and T' lines 
and at a generic u point for each phonon symmetry. For N 
even and N odd the selection rules are the same as for bi- and 
tri-layer graphene, respectively. 





BZ point 


phonon 


allowed 


scatt 


ering 


monolayer 


T(T') 


Ti 


Tb - 




T4 - 


T4 






T3 


T2 - 


-^n 








u 


u+ 


u~ 


— > u~ 






gated 


T(T') 


Ti 


Ti - 


-^Ti, 


Ta - 


T2 


monolayer 




T2 


Ti - 


^T2 








u 


u 


U — I 


u 






bilayer 


T(T') 


Ti 


Ti - 


-Ti, 


T2 - 


T2 


(N-cvcn) 




T2 


Ti - 


-^Ta 








u 


u 


u 


u 






biased bilayer 


T(T') 


T 


T - 


» T 








u 


u 


U — 1 


u 






trilayer 


T(T') 


T+ 


T+ 




T- 


T- 


(N-odd) 




T- 


T+ 


T- 








u 


u+ 


u+ 












u" 


u+ 


— > u~ 







VI. DOUBLE RESONANCE RAMAN 
SCATTERING PROCESS 

One example of explicit use of the electron-radiation 
and EPS selection rules is the double resonance Raman 
scattering process [s^ [s^, in which an electron in the 
conduction band is scattered by a phonon with wavevec- 
tor outside the F point in an intervalley (connecting elec- 
tronic states near the K and K' points) or in an intraval- 
ley (connecting electronic state near the same K or K' 
point) process. The G' Raman band (~ 2700 cm~^) 
comes from an intervalley process in which the electron is 
scattered by an in-plane transversal optic (iTO) phonon. 
We will discuss in details the G' scattering for mono- and 
multi-layer graphene. 

For the monolayer graphene, the possible scattering is 
illustrate in Fig. [3 The iTO phonon at the KM (T') 
direction presents a Ti symmetry 37[, which can only 
connect two electrons with the same symmetry. Many 
other similar scattering events are allowed by symmetry, 
involving electron in the KF (T) direction or at any gen- 
eral u point inside the circle defined by the T3 photon 
energy. However, the matrix element has a strong angu- 
lar dependence and the scattering is dominated by the 
T electrons, as discussed in Ref. [9]. Therefore, the G' 
Raman band has only one peak, with full width at half 
maximum (FWHM) of ~ 24 cmT^ (see Fig. [6] (a)) [l,!^. 
For the graphene on top of a substrate, the same selec- 
tion rules apply, and the expected number of G' peaks is 
the same as for the isolated monolayer graphene on an 



FIG. 4: The most efficient DRR process in graphene showing 
the absorption of hght with polarization symmetry T3 fol- 
lowed by electron scattering by a phonon with Ti symmetry. 




FIG. 5: The four allowed DRR processes in bilayer graphene, 
taking into account the optical anisotropy. 



isotropic medium. 

For the bilayer graphene, the number of allowed DRR 
processes predicted by group theory will be larger, since 
both electronic and phonon branches are doubled. Along 
the T line, there is polarization dependence for the ab- 
sorption linking different electronic bands, as discussed 
in section IIIIl For wavevectors in the range of visible 
light energy, the electron dispersion are almost linear, 
then optical anisotropy can be applied here as for mono- 
layer graphene (ssj . Now, for computing the number of 
resonant conditions involved in the DRR process, we are 
left with only two excited electronic bands with symme- 
tries Ti and T2, which corresponds to Fig. [3] (a). The 
iTO phonons for bilayer graphene have Ti and T2 sym- 
metries. For the electron scattering by a Ti phonon, the 
allowed process are between K and K' electronic bands 
with same symmetry (Ti — > Ti or T2 ^ T2). The same 
happens with the electron scattering by a T2 phonon, but 
it connects conduction bands of different symmetries, i.e. 
Ti T2. This gives rise to four possible DRR processes, 
as shown in Fig [5] Q] . The Raman spectra can then be 
used to differentiate mono- and bi-layer graphene (see 

Fig.imiii. 

In the case of biased bilayer graphene, there are no 
selection rules involving different photon polarizations. 



2550 2600 2650 2700 2750 2800 2850 

Raman shift (cm'^) 

FIG. 6: The measured Raman spectra of the G' band of a 
(a) mono- (b) bi- and (c) trilayer graphene. The samples 
were made by exfoliating graphite on top of a 100 nm silicon 
oxide substrate using 2.41 eV laser energy. The G' band for 
mono- bi- and tri-layer graphene were fitted with 1, 4 and 6 
Lorentzians, respectively, with a FWHM of 24 cm~^. 



The biased bilayer can have photon absorption linking 
all valence and conduction bands. This leads to eight 
possible transitions considering EPS selection rules for 
Ti and T2 phonons. 

For the trilayer graphene, the DRR process will have 
again more contributions because each phonon and elec- 
tron band will be split in three levels. Along the T line, 
there are five possibilities linking the electronic bands be- 
tween the K and K' points with a T"*" phonon, and four 
possibilities for the T~ phonon. The total number of 
DRR process predicted by group theory will be fifteen. 
However, the FWHM is large when compared to the en- 
ergy splitting between the G' Raman peaks, and when 
one makes measurements of the G' Raman band, these 
fifteen peaks cannot be distinguished, as illustrated in 
Fig. [51 Similar problem should happen for N>4. 



VII. SUMMARY 

In this work we analyzed the symmetry aspects re- 
lated to electrons and phonons at each point in the BZ 
of graphene, depending on the number of layers. The 
symmetry aspects can be generalized to any value of N, 
differing for N even or odd. For monolayer and bilayer 
we consider both an isotropic and an anisotropic medium. 
We derived the selection rules for electron-radiation and 
electron-phonon interactions. Some specific findings can 
be remarked: 

• For the monolayer graphene, the predicted optical 
anisotropy [28] comes out directly from group theory 
analysis. The electron-phonon scattering process is al- 
lowed by symmetry at any generic point (u) in the Dirac 
cone, and the observation of a single Lorentzian in the 
G' Raman band comes from a strong anisotropy in the 
electron-phonon matrix element Q. 
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• The gated graphene have lower symmetry, but the 
optical anisotropy is still present, and for the DRR 
process, the symmetry considerations are the same as 
graphene on an isotropic medium. 

• In the case of bilayer graphene, the optical anisotropy 
is also present and there are four dominant processes in 
the DRR. This number increases to eight on biased bi- 
layer. 

• In trilayer graphene, the number of possible DRR 
processes is fifteen. However, the 15 processes are not 
distinguishable and the G' Raman band can be nicely 
fit with 6 Lorentzians. Similar situation is expected for 
larger number of layers. 
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APPENDIX: NOTATION CONVERSION FROM 
SPACE GROUP TO POINT GROUP 
IRREDUCIBLE REPRESENTATIONS 

In this work we derived the and Tiat.vib for all points 
in the first BZ of multilayer graphene maintaining the no- 
tation of space group (SG) for the irreducible representa- 
tions. The conversion to point group (PC) representation 
is obtained considering that (a) superscript sign or 
"- applies if the character of the horizontal mirror (ah) 
or inversion (i) is positive or negative, respectively; (b) 
the subscript number is given following the order of the 
point group irreducible representations; (c) two repre- 
sentations can only have the same number if they have 
superscript with positive or negative signs. As an exam- 
ple we give in Table IVlTl the F point space group notation 
conversion to the D3/1 (N-odd) and D^d (N-even) point 
groups and for the K point space group to the C^h (N- 
odd) and D3 (N-even) point groups. 
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